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Abstract 

We show how to coarse grain polymers in a good solvent as single particles, interact- 
ing with density-independent or density-dependent interactions. These interactions 
can be between the centres of mass, the mid-points or end-points of the polymers. We 
also show how to extend these methods to polymers in poor solvents and mixtures 
of polymers. Treating polymers as soft colloids can greatly speed up the simulation 
of complex many-polymer systems, including polymer-colloid mixtures. 
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1 Introduction 



Binary mixtures of colloidal particles and non-adsorbing polymers have re- 
ceived renewed and growing attention recently, in part because they exhibit 
complex and interesting structure, phase behaviour, interfacial properties, and 
rheology [1-6], and in part because they are excellent model systems for the 
study of large length and time-scale separations in complex fluids. Problems 
with bridging length-scales are immediately apparent in even the simplest 
models of colloid-polymer mixtures: while the mesoscopic colloidal particles 
can be modeled as hard convex bodies, the polymers are generally treated at 
the microscopic (Kuhn) segment level. Thus, even though the average size of 
the polymer coils may be of the same order of magnitude as that of the colloids, 
the number of degrees of freedom needed to model the former may be several 
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orders of magnitude larger than what is needed for the latter. This naturally 
provokes the question: Can the polymers also be modeled as single particles? 
In fact, this is exactly what was done by Asakura and Oosawa (AO) who, in 
their classic work on colloid polymer mixtures[7], modeled the polymers as 
ideal particles with respect to each other, and as hard-spheres with respect 
to the colloids. This model is strictly speaking only valid for non-interacting 
polymers, or for interacting polymers in the dilute limit, while many interest- 
ing phenomena, such as polymer induced phase separation, take place at finite 
concentration of interacting polymers. Our ultimate goal, therefore, is to go 
well beyond the AO model and describe non-ideal polymers in a good solvent 
up to semi-dilute concentrations. We recently extended the AO concept to take 
into account polymer-polymer interactions, first by rather naively assuming a 
Gaussian repulsion between polymers [8] to account for the penetrable nature 
of polymer coils, and then by carrying out a more sophisticated programme 
which resulted in density-dependent [9-11] and density-independent [12] inter- 
actions between polymer coils. 

The next sections will examine these effective potentials in more detail. 



2 Coarse-graining homogeneous polymer solutions 

Polymers made up L segments are characterized by their radius of gyration, 
R g ~ L v , where v xs 0.59, i.e. polymers are fractal objects. Much of our 
understanding of polymer solutions comes from the scaling arguments pio- 
neered by de Gennes[13]. These arguments suggest that the behaviour of a 
polymer solution differs in the dilute regime, where the polymer coil density 
p < p* = §7T-Rg, from the semi-dilute regime, where p/p* » 1. Both these 
regimes presume that the actual monomer concentration c remains very low. 
Once c becomes appreciable, one enters the so-called melt regime which will 
not be treated here, and for which different kinds of coarse-graining methods 
are necessary. 

A standard scaling argument suggest that the radius of gyration R g is the 
only relevant lengthscale for the dilute and semi-dilute regimes of polymers in 
a good solvent [13]. This immediately implies that the second- virial coefficient 
should scale as B2 ~ Rg. If we set x = r/R g then the second- virial coefficient 
is proportional to 



where f3V(x) is the interaction between two separate polymers, defined w.r.t. 
some yet to be specified coordinate. Since this must hold for all R g , the in- 
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teraction f3V(r/R g ) should not depend on the length L for sufficiently long L 
(scaling limit). A more sophisticated version of this argument was put forward 
by Grosberg, Khalatur, and Khokhlov[14], see also the review by Likos[6] for 
a historical overview. 

To coarse-grain each polymer as a single entity, one must still choose an in- 
teraction centre, which may be the centre of mass (CM), the mid-point, the 
end-points, or some average monomer. For the mid-point or end-point rep- 
resentation, f3V(x) should diverge at the origin, since the actual segments of 
two different polymers cannot overlap. For the CM, we expect a finite value 
of /3V(x = 0), since it is possible for two polymers to deform around each 
other in such a way that their CM coincide without any mutually avoiding 
monomers overlapping. 

2. 1 Density-independent polymer-polymer interactions 

In the description of atomic and molecular liquids and solids it is common to 
replace the full quantum mechanical treatment of the interactions by a simpli- 
fied effective potential. Well known examples include the Lennard Jones pair 
potential and the Axilrod- Teller three-body potential [15]. Here we attempt 
a similar coarse-graining for polymer solutions and choose the constituents 
to be single polymers, with a specified interaction centre for each polymer. 
Then, following for example [5,6] or more specifically [11], the coarse-grained 
Helmholtz free energy F of a set of N polymers with their interaction centres 
fixed at the coordinates {r^}, in a volume V, can be written as the following 
expansion: 

N 

F(N,V,{r t }) = F^\N,V)+ £ ^ (2) (r n , r, 2 , ) (2) 

h<t2 

N 

+ w i3) {r h ,r i2 ,r i3 ) + . . . +w (N) {r h ,r i2 . . .r iN ) 

il<*2<*3 

where the coordinates of the interaction centres of the polymers, {r^, r i2 . . . r in }, 
are expressed in units of R g , the radius of gyration at zero density. Each term 
in this coarse-grained free energy includes an implicit statistical average over 
all the internal monomeric degrees of freedom for a fixed configuration {r^}. 
F(°)(iV, V) is the so-called volume term[6], the contribution to the free en- 
ergy that is independent of the configuration {rj}, and includes the internal 
free-energy of an isolated polymer. For a homogeneous solution, translational 
invariance implies that there is no one-body term in the expansion. Each sub- 
sequent term ■u/^r^, r i2 . . .r in ) is defined as the free energy of n polymers 
with their interaction positions at {r^, r i2 . . . r in }, minus the contributions of 



3 



all lower order terms. This procedure may in principle be followed to derive 
higher and higher order interactions, until, for a system with N polymers, the 
iVth term determines the total coarse-grained free energy. The thermodynamic 
free energy of the polymer solution finally follows from a statistical average 
over the interaction coordinates: 

PF(N, V) = - In ]T exp [-/3F(7V, V , {r,})] . (3) 



But in practice, this approach is not often feasible because the number of n- 
tuple coordinates and related complexity of each higher order term increases 
rapidly with n, so that the series in Eq. (2) and the full average in Eq. (3) 
quickly become intractable. Instead, one hopes to show that the series con- 
verges fast enough that only a few low order terms are needed to obtain a 
desired accuracy. 

The first important term in the series expansion is the pair interaction w^ 2 \r), 
which can be determined by calculating the logarithm of the probability that 
two polymers have their interaction centres a distance r apart. Details of our 
computer simulation technique are described elsewhere [10, 11, 16]. In brief, by 
simulating L = 500 self-avoiding walk (SAW) polymers on a cubic lattice, we 
determined w^(r) for three different interaction centres: the end-points, mid- 
point, and CM of each polymer, as depicted in Fig 2. The end and mid-point 
representations diverge at the origin as lim r ^ w^(r) ~ \n(r/ R g )[6], while 
the CM representation has a finite value which we estimate to be w^(0) = 
1.80 ± 0.05 in the scaling limit L — > oo[10]. By plotting r 2 v(r) we see that the 
CM representation has the shortest range, which is one reason why it is easier 
to use than the other two representations. 

In a similar fashion, the higher order interactions can be calculated from higher 
order probability distributions [5, 11]. We calculated the relative strength of the 
many-body terms up to fifth order by computer simulations, and to arbitrary 
order by a scaling theory [11]. The simulations and the scaling theory agree 
well, and suggest that at full overlap the nth order many-body term alternates 
in sign as (—1)™, and decreases (slowly) in absolute magnitude with increasing 
n. But, as mentioned before, a description based on three and higher order 
interactions rapidly becomes too unwieldy to use. Instead, we show in the next 
section how to derive a pair potential approach which includes these higher 
order (n > 2) interactions in an average way. 
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2.2 Density- dependent polymer-polymer interactions 

An alternative coarse-graining approach to the many-body expansion of the 
previous section is to find pair potentials which reproduce known structural 
information. We are aided in this by a theorem which states that at a given 
density p, there is a one-to-one mapping between the pair distribution function 
g{r) and a unique pair potential v(r; p) that will exactly reproduce the correct 
pair correlations [17]. We generate the pair correlations with simulations of 
SAW polymers, and at each density p, use the Ornstein-Zernike equations[15], 
coupled with the hypernetted-chain closure (HNC)[15], to invert the CM g(r) 
and find v(r; p). While the HNC closure is generally not accurate enough for 
inversions in simple liquids, it is nearly exact for the soft potentials we are in- 
vestigating here[5,6,18]. Nevertheless, there are a number of subtleties, both in 
the simulations and in the inversions, which must be carefully examined[10,12]. 

Density-dependent effective potentials v(r; p), inverted from the g(r) produced 
by L = 500 SAW simulations are shown in Fig. (3). The potential changes 
with increasing density, but approximately retains the shape found at p = 0. 
We have recently shown that these potentials can be very accurately param- 
eterized for p/p* < 2 by sums of three Gaussians with density-dependent 
coefficients [12]. 

Within the HNC approximation, the density dependence of an effective pair 
potential that reproduces the true g(r) is given to lowest order in p and the 
wW({Ti}) by[19]: 

v(r 12 ; p) = w {2 \r l2 ) -pj (V^m,™) _ ^ p )g 2 ^ p)dr 3 ,(4) 

We found that this expression describes the density dependence quite well 
for p/p* < 1, and even works qualitatively for higher densities, where we 
expect higher order p and w^ n \{ri}) effects to become significant [11]. This 
demonstrates the connection between the density-independent and density- 
dependent approaches, showing explicitly that the density dependence in the 
effective pair-potentials v(r; p) arises from the many-body interactions. 

One advantage of the structure-based route to the potentials is that one can 
use the compressibility equation[15] to derive the equation of state (EOS) 
from the pair correlations. We have done this for both L = 500 and L = 
2000 SAW simulations. We directly measured the EOS and compared this 
to the EOS derived from the effective potentials through the compressibility 
equation [10, 12]. The two routes are compared in Fig. 4, where the agreement 
is shown to be excellent. 

All three effective potentials shown in Fig 2 result in "mean-field fluids" [5,6,18], 
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so named because the EOS takes on the mean-field form /3H/p ~ l+pv(k = 0) 
at high enough densities. Here v{k) is the Fourier transform (FT) of the po- 
tential. This implies that if we only use the p = potential, then the EOS 
at higher density would scale as (3TL/ p ~ p instead of the the correct p ~ p L3 
scaling found for the semi-dilute regime[13]. It is therefore the many-body in- 
teractions, expressed through the density dependence of v(r;p), which cause 
the EOS to be super-linear. 

We add a caveat here about the route to thermodynamics with density depen- 
dent potentials. The potentials derived here can be used to derive the correct 
thermodynamics through the compressibility route. Different (but related) 
density-dependent potentials, which do not reproduce the correct structure, 
would be needed to derive the correct thermodynamics through the virial 
route [19]. In other words, there is no unique density-dependent pair poten- 
tial: when specifying such a potential, one must also specify which route to 
thermodynamics should be used [20]. 



3 Polymers near walls and spheres 



Polymers form a depletion layer near a hard non-adsorbing wall because the 
number of possible conformations are restricted there. This is illustrated in 
Fig 6 for the CM and for a monomer representation near a planar wall. We 
have used an inversion method similar to that used in the previous sections 
to derive effective wall-polymer potentials for interacting polymers near walls 
and spheres [9, 10, 12]. These potentials are constrained to give the correct den- 
sity profile p(r) which in turn determines the adsorption T. One can show, for 
example, that the surface tension is completely determined if one knows the 
EOS and T as a function of p[21]. Our effective polymer-polymer pair poten- 
tials correctly determine the EOS, while the wall-polymer potentials correctly 
determine T, implying that our formulation will reproduce the correct surface 
thermodynamics. 

We show this adsorption T/p in Fig. 7, together with a simple fit constrained 
to give the correct scaling T/p £(p) ~ p~°- 77 in the semi-dilute regime. Note 
that the largest relative change in the adsorption is actually in the dilute 
regime, suggesting that even there descriptions based on the low-density or 
non-interacting polymer limit rapidly become inadequate. 
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4 Connection with scaling theory 

Most successful theories of polymers start from a monomer based description 
and use scaling or RG approaches to derive properties of polymer solutions [13]. 
How does our CM based description compare with these scaling approaches? 
For example, in the semi-dilute regime, scaling theories predict that the im- 
portant length-scale is the correlation length £(p), which decreases with in- 
creasing density as £(p) ~ p~ a77 . It is not a-priori clear how this lengthscale 
enters into the g{r) or the v(r) in our description of homogeneous polymer 
solutions. The EOS scales as /311/p ~ £~ 3 /p in the semi-dilute regime, a be- 
haviour which is reproduced by our description through the compressibility 
equation. Since these potentials result in "mean-field fluids", this suggests 
that / drr 2 v(r; p) ~ £~ 3 /p 2 in the semi-dilute regime. For inhomogeneous sys- 
tems in the semi-dilute regime, £ enters more directly through the density 
profiles shown in Fig 6, but again the direct connection to the potentials is 
more opaque. So the exact connection with the scaling theory still remains 
to be worked out. We expect our approach to be most robust in the dilute 
regime and into the crossover region of the semi-dilute regime. Luckily this 
is also where much of the interesting physics of the colloid-polymer systems 
lies. How well our "soft colloids" approach will work deep into the semi-dilute 
regime still remains to be established. 



5 Extensions to poor solvents and mixtures 

The considerations in the previous sections focused on equal length polymers 
in a good solvent, where the temperature plays no role. However, the tech- 
niques used for polymers in a good solvent should still apply to other types of 
polymers solutions. 



5. 1 Poor solvents 

We first examine briefly what happens for polymers in a poor solvent, using as 
a model SAW polymers with a nearest neighbor attraction of strength —f3e. 
It is known that as the temperature decreases, there is a temperature T co i 
below which the polymer collapses into a compact globule and loses its fractal 
nature [13]. The effective potentials will then be fundamentally different of 
course. But as long as we stay above this temperature, we expect that the 
interaction should become less strong with decreasing temperature, as shown 
in Fig. 5. 
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5. 2 Mixtures 



Renormalization group (RG) calculations for the interaction between the CM 
of two polymers of differing lengths R g i and R g 2 suggest that the interac- 
tion strength at full overlap should weaken with increasing size asymme- 
try, and that the interaction range should approximately scale as R12 = 



\\J R 2 g i + -R32P2]. We confirm this behaviour for simulations of a number of 
different length SAW polymers at p = in Fig. 5. 

5.3 Phase separation? 

An interesting prediction is that binary mixtures of Gaussian core particles 
interacting via the pair potentials V a/ 3(r) = e aj gexp{ — (r/R a p) 2 } will phase- 
separate over broad ranges of the coupling constant ratio 1 612 1 / y/euOw ■ Semi- 
quantitative correct phase-diagrams follow already from an analytic mean-field 
calculation [23], and suggest that polymers in a good solvent will not phase 
separate at low densities. 

5.4 Relationship with PRISM 

Each time the present coarse-graining methods are applied to a new type of 
polymer solution we need a new set of computer simulations at the monomer 
level for the parameters of interest. It would be very helpful to find other, 
semi-analytic, ways of providing input information that are faster and more 
flexible. One candidate would be PRISM [24], an integral equation method 
which has been applied to a wide variety of polymeric systems. This requires 
a way of deriving CM-CM correlations from the monomer based correlations 
provided by PRISM. As a first step in this direction, we have derived an 
approximate relationship which is much more accurate than earlier, heuristic 
approaches [25]. This could form the basis for using PRISM or other monomer 
based methods as input to our "polymers as soft colloids" approach. 



6 Conclusions 

In summary then, we have shown how to coarse-grain polymers as single "soft 
colloids", with just three degrees of freedom each, interacting via a density- 
independent pair, triplet, and higher body potentials. These, however, become 
rapidly intractable. In parallel, we also derived density-dependent pair poten- 
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tials which include, in an average way, the effect of the higher n-body inter- 
actions. These effective pair potentials exactly reproduce the two-body corre- 
lations of the underlying polymer solution, and, through the compressibility 
equation, they reproduce the EOS as well. In a similar way the effective one- 
body potentials which exactly reproduce the one-body density profiles near 
walls and spheres can also be derived. Because these reproduce the correct 
adsorptions, the thermodynamics of a polymer solution near a non-adsorbing 
surfaces are also correctly reproduced by our formulation. 

We also showed how to extend these methods to derive effective potentials for 
polymers in a poor solvent and for mixtures of different length polymers. We 
sketched some ways in which other monomer based methods such as PRISM 
could be used as the source of input to derive our potentials. 

At this point one might ask what has been gained, since at each point di- 
rect computer simulations were needed as input to derive the potentials. This 
question brings us back to the aim stated at the outset: to describe mixtures 
of many polymers and many colloids. Here our coarse-graining of polymers 
as soft colloids does result in important simplifications. For example, we have 
performed such simulations for mixtures of spheres of radius R c and poly- 
mers with sizes ranging from R g /R c 0.3 to R g /R c ~ 1, and determined the 
polymer induced phase-separation of the colloids[26]. The effective polymer 
density p/p* at the critical point increases with increasing polymer size, but 
even for the largest polymers it is still in the regime p/p* < 1, where we expect 
our formulation to work best. Simulations with a full polymer model would 
be about two orders of magnitude slower. Without coarse-graining the poly- 
mers as soft colloids, such a simulation would have been virtually impossible 
to perform. 
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Fig. 2. w( 2 \r) and r 2 w^ 2 \r) for the interaction between two isolated polymers in 
the CM, end-point, and mid-point representations. The end-point and mid-point 
potentials diverge at the origin, but the CM representation gives a finite value. All 
three potentials result in the same second virial coefficient B2 defined in Eq. (1). 




Fig. 3. The effective polymer pair potentials v(r;p), derived at different densities 
from an HNC inversion of the CM pair distribution functions g(r) of L = 500 SAW 
polymer coils, (from Ref[ll]) 
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Fig. 4. EOS for polymers in good solvent. Direct results and the EOS arising from 
the effective potentials through the compressibility route are compared for L = 500 
and L = 2000 SAW polymers. 
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Fig. 5. (a) CM-CM interaction for polymers in a poor solvent, modeled by L = 100 
SAW polymers on a cubic lattice with nearest neighbour attractions of strength 
—fie. (b) CM-CM interaction for mixtures of different length SAW polymers. 
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Fig. 6. (a) The wall-polymer CM profile h{z) = p{z)/p - 1 for L = 500 SAW 
polymers at different bulk concentrations, (b) The wall-polymer monomer profile 
h m (z) for the same bulk concentrations. Both representations result, by definition, 
in the same relative adsorptions T/p. The straight lines in (b) are a fit to the simple 
form h m (z) = tanh 2 (-zp/T(p)) - 1[13] 
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Fig. 7. Relative adsorption T/p, in units of R g of L = 500 SAW 
polymers near a single hard wall. Circles are direct simulations and the 
line denotes the simple fit with the correct scaling behavior, given by 

/ „\ -(0.2565) 

T{p) = -l.WARg (l + 7.63^ + 14.56(^) 3 J 
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